Segregation pattern reorientation of granular mixture on horizontally oscillating tray 

Masashi Fujiij*] Akinori Awazu, and Hiraku Nishimori 

Department of Mathematical and Life Sciences, 
Hiroshima University, Kagamiyama, 
Higashi- Hiroshima, 739-8526, Japan. 



(N 

o 

(N 

< 

in 

CM 

d 



Reorientation of the segregation pattern of a binary granular mixture on a two-dimensional hor- 
izontally oscillating tray is numerically realized. The mixture consists of large-and-heavy particles 
and small- and-light particles, the segregation pattern of which shows a transition between a striped 
pattern perpendicular to the oscillation and that parallel to the oscillating direction according to 
the change of area-fractions of two types of particles. The transition mechanism is discussed on the 
basis of a simplified 1-dimensional dynamics. 
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Granular materials exhibit various complex behaviors, and their characteristic dynamics have been studied exten- 
sively, for instance, the crystallization of granular materials l| or the pattern formation of granular mixtures on a 



horizontally oscillating tray 



0Q. 



2|-l5( , in a horizontally rotating drum 6Mll| , and in a vertically oscillating box 
However, thus far, the research on the basic physics of the granular dynamics as well as macroscopic phenomenology 
has been insufficient. For example, previous studies on the patten formation of a granular mixture on a horizontally 
oscillating tray 2|-|5| have demonstrated that different types of granular materials segregate to form a stripe pattern 
that extends perpendicular to the oscillating direction, while no segregated pattern has been observed parallel to the 
oscillating direction. The potential possibility of the emergence of the latter segregation pattern has been an intriguing 
issue in this field of study. For a wider range of soft materials, a recent numerical study on a colloidal mixture found 



14j. In this study, we use a 



stripe patterns that were both parallel and perpendicular to the oscillating electric field 
numerical model and show that transition between steady perpendicular and steady parallel stripe patterns can be 
observed in granular mixture as well, for a suitable set of parameters and a suitable combination of the area fractions 
of the two types of constituting particles. 
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We employ the distinct element method (DEM) 3|, |4| and simulate the behaviors of the granular mixture on a 
two-dimensional horizontally oscillating tray. This mixture consists of two types of particles: i) the particle with a 
large grain size and mass, which we refer to as the L-particle, and ii) the particle with a small grain size and mass, 



which we refer to as the S'-particle (see Table H]). The i-th particle has type-dependent mass (i 



types 



and size (D t i ype ) 



(type = L or S), and its position is denoted by a two-dimensional vector, r,(= (xi,yi)). Each particle is subjected 
to viscous friction, which is proportional to the velocity difference between the particle and the tray with the friction 
coefficient [i 3|, If the i-th and j-th particles overlap, the interaction is indicated by a linear elastic force with 
viscous dissipation, depending on the combination of the types of particles; the elastic coefficient and the coefficient 
of dissipation are denoted as k and 7^, respectively. Here, 7^ is derived from the coefficient of restitution e l ' J , the 
elastic coefficient k, and the reduced mass of particles as follows. 
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The dynamic equation of the i-th particle is summarized as follows. 
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(for \6r\ < (Dl vpe + Df pe )/2) 
(for \St\ > {D* ype + Df pe )/2) 



(1) 



(2) 



(3) 



where the velocity of tray is given as 



5r itj = r, - r 
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b tmy = (*tray>ytray) = ( 27rAu sm(2-Kut) , 0) 
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Note that the rotation of the particle is not considered. The values of parameters are summarized in Table HI most 
of which are set same or close to those introduced in previous related studies 3|, As one of exceptions, the elastic 
coefficient k is set smaller than that in previous study, meaning that the particles are not so rigid, but the collision 
time derived from the elastic coefficient remains small compared with the period of oscillation of tray. In addition, for 
the simpleness, we consider only the case where the elastic coefficients of L- and S'-particlcs arc identical. Another 
exception is the large value of /1 assuming that particles are more strongly affected by the viscous friction arising 
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TABLE I: Parameter values considered in numerical simulation for pattern transformation. 



property 


L- value 


S-value 


m type : Mass 


m L = f g 


m s = 0.025 g 


D type : Size 


D L = f cm 


D s = 0.5 cm 


(L,S) or (S,L) 
(S,S) 


e L ' L = 0.2 
e L ' s = 0.5 
e s ' s = 0.9 


H : Friction coefficient 


f gs- 1 


k : Elastic coefficient 


f.O x fO 4 g cm 2 s -2 


A : Amplitude of oscillation of tray 


fO cm 


v : Frequency of oscillation of tray 


2 Hz 



from their relative velocity to the tray than previous studies. This is to explore the potentially new aspect of the 
segregation dynamics on a horizontally oscillating tray. Here, the random fluctuation in the system is neglected. 
Systematic calculations to verify the effect of viscous friction on segregation will be presented later. 

We perform simulations using several combinations of Nl and Ns, «-e. the number of L- and S-particles, respectively. 
In addition, wc define the area fractions of L- and S'-particles as pl = Nltt(D l ) 2 / {AW 2 ) and ps = Nsn(D s ) 2 /(AW 2 ), 
respectively. The width of the system is W = 40 cm in both oscillational and lateral directions, and the periodic 
boundary conditions are set in both directions. Each run starts from t = 0, when the particles are randomly distributed 
on a tray without overlap, and ends at t = 200. 

Figure [1] shows typical patterns emerging as steady states of the binary granular particle mixtures: (a) non- 
stripe pattern, (b) stripe pattern perpendicular to the oscillating direction (hereafter, we refer to this pattern as 
the perpendicular-striped pattern), (c) stripe pattern parallel to the oscillating direction (hereafter, we refer to this 
pattern as the parallel-striped pattern), and (d) crossed-striped pattern. 

To distinguish these four patterns quantitively, we define a temporal order parameter $(i). Specifically, the system 
is separated into M x M (M = 40) square cells and a quantity, 

= N&t) ■ (D L ) 2 - Nfjt) ■ (D s ) 2 , (6) 

is allocated to each cell (i,j), where N^(t) and Nfj(t) are, respectively, the number of L- and S'-particles in cell (i,j) 
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FIG. 1: Typical steady states of particle mixtures: (a) non-stripe pattern (pL = 0.1, ps = 0.65). (b) stripe pattern perpen- 
dicular to oscillating direction (p_L = 0.30, ps = 0.55), (c) stripe pattern parallel to oscillating direction (pL = 0.30, ps = 0.25), 
and (d) crossed-striped pattern (pL ~ 0.40, ps = 0.40). The scale bar denotes 20cm. 

at t. Then, the specific form of the temporal order parameter $(f) is given as 



The first and second terms of the right-hand-side of Eq. ([7]) quantify the spatial correlations along the direction 
parallel and perpendicular to the oscillating direction, respectively. In simple terms, approaches 1 when the 
parallel-striped pattern is dominant, whereas it approaches -1 when the perpendicular-striped pattern is dominant. 

Simulations are systematically performed by varying the combination of ph E {0.05,0.1,0.15,- • -0.8} and ps € 
{0.05, 0.1, 0.15, • • -0.8}, maintaining ps + PL < 0.85. The relationship between (pl,Ps) and with averaged 
over the last 20 periods of oscillation, for A = 10 and v = 2 is shown in Fig. [5] As shown in the figure, to observe this 
relationship, wc consider four typical regions: (a) p^ < 0.1 or ps < 0.1, where $ a „ e is close to zero; (b) pl+Ps < 0.85, 
p L > 0.15, and p s > 0.45, where $ au is close to -1; (c) pl + Ps < 0.55, 0.15 < p L < 0.40, and 0.15 < ps < 0.35, 
where <E> at) is close to 1; and (d) 0.55 < pl + Ps < 0.85 and 0.20 < ps < 0.40, where <& al , e is close to zero. In 
region (a), particles tend not to form any stripe pattern, and in regions (b), (c), and (d), particles tend to form the 
perpendicular-striped, the parallel-striped, and the crossed-striped patterns, respectively. 





(7) 
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FIG. 2: Relationship between the order parameter <J> a „ and the area fraction of i-particles (horizontal axis) and S'-particles 
(vertical axis), for A — 10 and v = 2. A bright color implies that $„ is close to 1, indicating the emergence of a parallel-striped 
pattern, and a dark color implies that $ a „ is close to -1, indicating the emergence of a perpendicular-striped pattern. 




Hi 



FIG. 3: Relationship between the order parameter & av and the friction coefficients of L-particles (horizontal axis, log scale) 
and S'-particles (vertical axis, log scale), for pL = 0.25, ps = 0.3, A — 10 and v = 2. A bright color implies that <& av is close 
to 1, indicating the emergence of a parallel-striped pattern, and a dark color implies that <& a v is close to -1, indicating the 
emergence of a perpendicular-striped pattern. 



It should be noted that the parallel-striped pattern obtained in region (c) has not been realized in previous studies 
on the horizontally oscillating granular mixture, wherein much lower viscous friction (/x ~ 0.1) has been used. To 
verify the effect of high viscosity friction, pi = 1, considered in the above calculations, as shown in Fig. [3j we calculate 
the order parameter in steady state with a varying the combination of friction coefficient of L-particles and S'-particles 
keeping pl = 0.25, ps = 0.3, A — 10 and v = 2. The case of fiL = fJ-s = 1 in this figure corresponds to the case 
of pl = 0.25, ps = 0.3 in Fig. [21 As roughly shown in Fig. [31 particles tend to form the parallel-striped pattern in 
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FIG. 4: (a)Time evolution of the order parameter, Q(t), for pl = 0.30 and ps = 0.25, and the corresponding simulation 
snapshots at (b) t = 0, (c) 30, (d) 60 and (e) 70. The scale bar denotes 20cm for (b)-(d). 

the high-viscous-friction range, and the perpendicular-striped pattern in the low-viscous-friction range. This is the 
reason why we choose the set, hl — Ms — 1> °f viscous friction coefficients in this study with which set we can see 
the present characteristic segregation dynamics approaching the parallel-striped pattern. Now, we focus on the time 
evolution of $(f) in region (c). Figure BJb)-|3|e) shows simulation snapshots at t = 0,30,60, and 70, respectively, 
at (pLiPs) = (0.3,0.25); there, the gray particles correspond to L-particles and the black particles correspond to 
5-particles. At the early stage of the simulation, < t < 30, particles temporarily form the perpendicular-striped-likc 
pattern as shown in Fig. where the value of $(t) is kept around zero. In this stage, the perpendicular-striped 

pattern arises, but alignments of stripes are not so distinct. In the next stage t = 30 until t = 90, the value of <&(£) 



%f%2 "' X Ns X l X 2 "' X N L 
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FIG. 5: Illustration of the system of a one-dimensional model. The small black circle represents the S-particle and the large 
gray circle represents the L-particle. 

gradually increases, accordingly, the striped-pattern is largely deformed to form clusters of L-particles. For example, 
after t — 60, the perpendicular stripe is hardly recognizable as seen in Fig. EJd), instead, an wedge of horizontally 
extending cluster of S-particles is seen to penetrate into the the clusters of L-particles, which serves as a seed of parallel- 
striped pattern. Subsequently another horizontal wedge of S-particles begins to collide the clusters of L-particles to 
form the 2nd parallel stripe as seen in FigUJ^e) for t = 70. Around t = 90, $(i) reaches a saturation value &(t) ~ 0.55, 
before which time the parallel-striped pattern has been established, in that situation, the perpendicular movement 
of particles is restricted by the strong friction force acting between the particles and the horizontally oscillating tray. 
Hence, the movements of L- and S-particles get synchronized, resulting in a sharp decrease in the frequency of particle 
collisions. In this way, the deformation of perpendicular stripes, which are formed once in the system, irrespective of 
the final pattern, plays a key role in the selection of the final pattern. In particular, considering the far more active 
movement of S'-particles than L-particlcs, that causes the deformation of the perpendicular alignment of L-particlcs, 
the stability of stripes of L-particles against the collisions of the cluster of S-particles seems essential for the selection 
of final pattern selection; the perpendicular-striped pattern or the parallel-striped pattern. 

Accordingly, to investigate the stability of L-particlc stripes (cluster) in a simplified setup, we consider a one- 
dimensional (ID) granular system, in which particles arraying in the x (oscillating) direction share the central axis, as 
shown in Fig. [5j the L- and S'-particles are kept unmixed. This is the ID expression of a-pair-of-perpendicular-striped 
pattern realized in the 2D system. The dynamics and external force are kept the same as those in 2D system, where 
the boundary condition is periodic. Note that because of the dimensionality, the sequence of particles is preserved in 
the ID simulation. In this situation, we do not focus on the stability of the initial sequence of the L- and S-particles 
but on the stability of the agglutinated state of the L-particles. This agglutinated state of the L-particles corresponds 
to a rigid structure of individual perpendicular stripes. 

To quantify the agglutination degree of the L-particles in the steady state of the system, we measure the mean 
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FIG. 6: Relationships among Nl (horizontal axis), Ns (vertical axis) and ij„ e , and MSD. The very dark color indicates that 
the value of x ave is smaller than 1. The lighter color indicates that the value of i'„ e is larger than 1. 

distance between neighboring L-particlcs, 

l _ 1 [ T KM - *m modJWQ 

T J T N L -1 

where W(— 20) is the total length of the system. 

The relationship among Nl, Ns, and x^ ve , obtained by the simulation, is shown in Fig. [6] The very dark color 
indicates that the value of x ave is smaller than 1, the lighter color indicates that the value of x% ve is larger than 1. 
In this figure, we see three typical regions. In the cases where Ns < 15 and Nl < 9, the mean distance between the 
L-particlcs is larger than 1; in other words, the L-particles cannot be kept agglutinated in the field. It should be 
noted that, in this region, L- and S-particles frequently collide at their boundaries. In contrast, in the case where 
Ns > 20, the obtained values of x% ve are smaller than 1, implying that the agglutinated state of L- particles is stably 
maintained. In this case, the frequency of the collisions between the S- and i-particles is much lower than that 
in the above case. In the case where Nl > 10, x^ ve takes intermediate values close to 1. These results suggest 
that the destabilization of the agglutinated state of L-particles in ID simulation is realized for small spatial densities 
of L- and S'-particles; this destabilization of the agglutinated state roughly corresponds to the unstable condition 
of perpendicular stripes of L-particles in 2D simulation shown in Fig. [21 On the other hand, in the high-density 
region of the ID and 2D simulations, particularly in the high-density regime of S'-particles, the agglutinated state 
of the i-particles is maintained stable. In this way, the formation of the parallel-striped pattern is caused by the 
destabilization of the L-particle stripes. This destabilization requires frequent collisions of S-particles, that is, the 
momentum transfer of S-particles, to the L-particlcs striped structure. This momentum transfer is realized when 
S-particles move actively in the field for a proper combination of the area fractions of S- and L-particlcs. The two 



destabilization processes are to be investigated in greater detail in future research. 

We numerically investigated the pattern formation of granular mixtures on a horizontally oscillating tray and found 
that a stripe pattern parallel to the oscillating direction is formed for small area fractions of particles. Focusing on 
the formation process of the parallel-striped pattern, the particles temporarily form a stripe pattern perpendicular to 
the oscillating direction of the tray in the early stage, hence, the perpendicular stripes of large and heavy particles 
are broken by repeated collisions with actively moving small and light particles, finally resulting in the formation and 
stabilization of the parallel-striped pattern^ 

It is noticeable that Alam et al. (2008) 



151 ] reported the emergence of parallel-striped patterns in the horizontally 
oscillating particulate suspension. In the present study, we did not consider the fluid motion, however, these two 
systems potentially belong to the same universality class in terms of instability of the striped pattern formation 
parallel to horizontally oscillating external force. The common backgrounding mechanism for the emergence of this 
type of patterns should be explored as the next issue. Also the more systematic study for the dependency of the final 
segregation pattern on the elastic constants of L- and S'-particle remains to be made. 
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